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Abstract 

We study exact effective superpotentials of four-dimensional M = 2 supersym- 
metric gauge theories with gauge group U(N) and various amounts of fun- 
damental matter on K 3 x S 1 , broken to N = 1 by turning on a classical su- 
perpotential for the adjoint scalar. On general grounds these superpotentials 
can easily be constructed once we identify a suitable set of coordinates on the 
moduli space of the gauge theory. These coordinates have been conjectured to 
be the phase space variables of the classical integrable system which underlies 
the M = 2 gauge theory. The sought low energy effective superpotential can 
then be constructed from the conserved quantities in the integrable system. 
For the gauge theory under study these integrable systems are degenerations of 
the classical, inhomogeneous, periodic <SX(2,C) spin chain. Ambiguities in the 
degeneration provide multiple coordinate patches on the gauge theory moduli 
space. By studying the vacua of these superpotentials in several examples we 
find that the spin chain provides coordinate patches that parametrize holomor- 
phically the part of the gauge theory moduli space which is connected to the 
electric (as opposed to magnetic or baryonic) Higgs and Coulomb branch vacua. 
The baryonic branch root is on the edge of some coordinate patches. As a prod- 
uct of our analysis all maximally confining (non-baryonic) Seiberg-Witten curve 
factorizations for Nf < N c are obtained, explicit up to one constraint for equal 
mass flavors and up to two constraints for unequal mass flavors. Gauge theory 
addition and multiplication maps are shown to have a natural counterpart in 
this construction. Furthermore it is shown how to integrate in the meson fields 
in this formulation in order to obtain three and four dimensional Affleck-Dine- 
Seiberg-like superpotentials. 
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1 Introduction 



The determination of the precise structure of gauge theories in the infra-red is an 
important and elusive problem. Recently, much progress has been made in calcu- 
lating the exact superpotentials in the low energy effective action in a large class of 
minimally supersymmetric (TV = 1) four-dimensional gauge theories. These theories 
are obtained by breaking A/" = 2 to M = 1 by turning on a classical superpotential. 
One way to study these theories involves the relation between M = 2 theories and 
integrable systems: the Seiberg-Witten curve of M = 2 theories can be identified with 
the spectral curve of an underlying integrable system. This identification relates the 
gauge theory moduli and the conserved charges (action variables) of the integrable 
system. If the gauge theory is compactified on R 3 x the relationship between 
gauge theory and integrable system becomes tighter, since the extra moduli from the 
compactification can be identified with the angle variables of the integrable system. 

Seiberg and Witten pQ have shown that the moduli spaces of four dimensional 
M = 2 gauge theories on K 3 x S l have a distinguished complex structure which is 
independent of the radius of the S 1 . The superpotential, which is a holomorphic 
quantity, will therefore also be independent of the radius of the S 1 , and the same is 
true for the vacuum structure of the theory. They do however change in the three- 
dimensional limit R — > 0, since in that case we should work with the 3d gauge coupling 
instead of the 4d scale A, and the relation between the two explicitly involves R. In 
this paper we restrict attention to finite R, as the generalization to zero radius is 
rather straightforward. 

Though the vacuum structure of the gauge theory on M 4 is the same as that 
on M 3 x S 1 , the field content and in particular the superpotential of each are quite 
distinct. There is a significant advantage to working on M 3 x S 1 compared to R 4 . 
The four dimensional gauge field gives, after compactification on S l , rise to a pair of 
scalar fields in three dimensions. We can give these scalar fields a vacuum expectation 
value which allows us to study the theory effectively at weak coupling. In particular, 
all non-perturbative effects will be due to conventional three-dimensional instantons. 
This in contrast to the four- dimensional situation where the non-perturbative physics 
is due to more complicated gauge field configurations such as fractional instantons. 

Based on these considerations and on previous work [21 El El, a precise conjecture 
was made in [Hj about the sought exact superpotential on K 3 x S 1 : it can be obtained 
by replacing the gauge invariant operators appearing in the classical superpotential by 
corresponding conserved quantities of the underlying integrable system. To achieve 
this, one can simply replace the adjoint superfield appearing in the superpotential 
by the Lax matrix of the integrable system. This is a very suggestive operation 
reminiscent of a master field, though a precise interpretation along these lines has 
not yet been found. In [H] this conjecture was verified by calculating the extrema 
of the resulting superpotential in various examples, and a precise agreement with 
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known results in four dimensions was found, and in particular the Seiberg-Witten 
curve always factorized in the appropriate way. A general proof was given in [01 [Zj - 
Several generalizations have been made since, to M = 1* [7 , to other classical gauge 
groups [H] and even to theories with gauge group G 2 0- In [HI] SU(2) with four 
flavours was studied. 

In this paper we discuss the generalization to U(N) gauge theories with Nf hyper- 
multiplets in the fundamental representation of the gauge group. According to [llj, 
the integrable system for Nf < 2N C should be a degeneration of the classical version 
of the inhomogeneous XXX spin chain. It was shown in that paper that the form of 
the spectral curve (for Nf = 2N C ) matched that of the Seiberg-Witten curve as in 
|12j . Quite generally, any holomorphic integrable system with the property that (i) 
the conserved charges (the action variables) precisely parametrize the moduli space 
of the four-dimensional gauge theory, (ii) the spectral curve always agrees with the 
Seiberg-Witten curve, and (iii) the angle variables parametrize the Jacobian of the 
Seiberg-Witten curve, is a candidate integrable system that one can use to describe 
the quantum superpotential of the gauge theory on I 3 x In this case the canonical 
variables (e.g. the coordinates and momenta) of the integrable system are complex 
coordinates on the moduli space of the gauge theory on M 3 x S 1 , and the quantum 
superpotential is simply an appropriate linear combination of the conserved charges. 
Thus, the main purpose of the integrable system is to provide a suitable set of complex 
coordinates on the moduli space. Unfortunately, these coordinates are not generically 
good global coordinates. If they were, the equations that put the conserved charges 
equal to fixed numbers would define a torus (the Jacobian) embedded in the phase 
space of the integrable system. Except for the one-dimensional torus, which can be 
written as a cubic equation in C 2 , the equations that describe complex embeddings 
of complex tori are very complicated. Therefore, what we will find is that the inte- 
grable systems parametrize an open subset of the Jacobian but not the full moduli 
space of the gauge theory on R 3 x S 1 . Correspondingly, we will see that some vacua 
of the gauge theory are not captured by the integrable system, or that they live at 
the boundary of phase space. Curiously enough, the 'electric' vacua of the theory 
are typically all recovered, while the dual 'magnetic' vacua are not. Another phe- 
nomenon that we will encounter is that there are several different integrable systems 
all describing the gauge theory with Nf < 2N C . These arise because there are several 
different ways in which we can degenerate the XXX spin chain. The different inte- 
grable systems parametrize different open subsets of the Jacobian, and on the overlap 
the relation between them is given by a suitable canonical transformation between 
the phase space variables. 

This paper is structured as follows: In section [21 we discuss the gauge theory 
expectations for the vacuum structure of theories with matter in the fundamental 
representation. Section 01 proceeds by introducing the classical spin chain, its spec- 
tral curve and its degenerations (down to eventually Toda). The following section, [3J 
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describes the natural proposal for the exact superpotential. The addition and mul- 
tiplication maps are discussed in section EJ Section El contains the result of several 
example calculations and some words to clarify their meaning. Specifically, we show 
how particular degenerations can sometimes lead to extra vacua. In addition, here 
all massive (field theory) vacua for a generic superpotential are obtained. In section 
Qit is shown how to integrate in meson fields into our proposed superpotential. We 
end with a discussion and present some conclusions and open problems. In several 
appendices we present some technical details, including the resolution of a point of 
minor confusion in the literature for high order superpotentials and an observation. 

2 Gauge theory expectations 

In this section we will very briefly review the zoo of phases of (the low energy effective 
theory corresponding to) Af = 2 U(N C ) 4 dimensional gauge theories with Nf (< 
2N C ) fundamental matter hypermultiplets, softly broken to Af = 1 by turning on a 
superpotential (the literature on this subject is vast, see e.g. [TH], [Mj or the review 
in J3| and references therein). 

2.1 Supersymmetric gauge theories on M 4 

The perturbative field content of the Af = 2 supersymmetric gauge theory we will 
be studying consists of the U(N) gauge vector multiplet, which contains an Af = 1 
adjoint scalar superfield and Nf pairs of massive hypermultiplets (Q, Q) in the 
(anti) fundamental. The global classical symmetries are an 577(2) x U (1) R-symmetry 
and an SU(Nf) x SU(Nf) xU(1)b flavour symmetry in the limit of vanishing masses. 
The superpotential for theories with Af = 2 supersymmetry and fundamental matter 
is 2 

VtV =2 = + V2Q?m}Qi, (1) 

ml — Diag(mi,m 2 . . . ,m Nf ). (2) 

It is well known that this theory is perturbatively asymptotically free in the UV for 
Nf < 2N C and finite for Nf = 2N C . We will break Af = 2 softly to Af = 1 by turning on 
a polynomial superpotential W n +i(<p) = J2k=o ^tr($ fc ) for the Af = 1 adjoint scalar 
$. These theories have a rich vacuum structure, which can be studied by a variety of 
methods. Below we briefly present the results of a field theoretical analysis. The first 
step is calculating the minima of the classical gauge theory, in order to identify the 
possible phases of the theory. Since some of them will not receive quantum corrections 

: We will denote the scalar components of the multiplets by the same symbol 
2 With JSj we follow the conventions of 
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the results there are exact. The D-term equations of the theory read 



[$,$ f ] = 

QW) b i - QXQVi = <, ^i. (3) 

with v the Fayet-Illiopoulos parameter, which in almost all cases must be zero in 
order for the classical theory to have a supersymmetric vacuum. We will not consider 
non-zero v in this paper. The F-terms give 

W'(<Z>) b a + V2Q°Qi = (4) 
<S> b a Ql + m)Qi = (5) 
Q b ^t + Q» = (6) 

As is well known, the first D-term equation states that gauge transformations can be 
used to diagonalize $ = Diag(0i, 2 , 3 , . . - 4>n c )- The F-terms state that the quark 
fields are either zero or eigenvectors of <£> with <p a = —mi for some a and i. If the 
quark masses are all distinct, we can label the quark fields and the adjoint scalar 
eigenvalues for which this condition is satisfied by a — i — 1, ... ,r. For every r there 
are ( Nf ) ways of choosing masses for which this condition holds. For unequal masses 
there are then r distinct eigenvectors of $ which we can group into diagonal 3 N c x Nf 
matrices for Q and Q*, 

Q = Diag(Qi,...,Q r ,0,...) (7) 
£? = Diag(Qi,...,g r ,0,...). (8) 

If some masses happen to be equal, than instead of a diagonal Q matrix we get a block 
diagonal matrix Q with sub-blocks Qi, but we will first assume unequal masses. The 
remaining D and F term equations now split into equations for the different diagonal 
elements in Q. We get 

IQ,l 2 = l4l 2 . (9) 

Next, we order the masses and the eigenvalues of the adjoint scalar such that the first 
r of each are equal. The value of the quark vevs and the remaining eigenvalues of the 
adjoint scalar can be determined from the remaining F term equations, 

W'(-mi) + V^QiQ* = 0, i = 1, . . . , r (no sum) (10) 
W'(cf) a ) = 0, a = r + l,...,N c . (11) 

If W'(-mi) is zero for some i, the quark fields Qi should be zero and do not give 
rise to a Higgs branch. If W'(-mi) is non-zero, the quark gets a vev which Higgses 
the corresponding U(l) gauge group. Clearly, the number of vacua depends on the 

3 by diagonal we mean = if i ^ j and if i, j > min(iV c , Nf) 
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detailed form of the superpotential. For example, a quadratic superpotential will 
force all the N = N c — r eigenvalues of the adjoint scalar <fr a which are not equal to 
some —rrii to be equal to each other. The resulting U(N) theory is known to have N 
super symmetric vacua. 

When some masses are equal (as they will be in most of the cases we consider) 
more elaborate vacua are possible. See the analysis in, for instance, ^3]- In summary 
the following happens: If iVj eigenvalues of the adjoint become fixed at one mass we 
retain a U(Ni) gauge theory with a number Mj of massless flavors. If W'(— mi) is not 
zero, the combination of equations Q and (fTTj) (suitably generalized for matrices Q 
and Q) forces Q and Q to acquire iVj non-zero eigenvalues which completely Higgs 
the U(Ni) theory. If W'(—m,i) is zero, at least some Q and Q must be zero (and, in 
fact, they may all be). The solutions are then labeled by an index I < min(-^, Ni) 
and generically the quark fields both have I degrees of freedom, constrained by the 
remaining D term equation. The theory that described the remaining degrees of 
freedom is then a U(Ni — I) gauge theory with M, — 21 massless flavors and l[Mi — I) 
massless neutral Goldstone hypermultiplets from the broken flavor symmetry. 

Summarizing, the theory has several phases, which we distinguish by the expecta- 
tion values of the quark fields. If these are zero, then we are on the Coulomb branch. 
If they are nonzero then we are on the Higgs branch. It can happen that the expec- 
tation values of the quark fields on this branch can be tuned to zero, in this case we 
can continuously interpolate between Coulomb and Higgs branch. The special point 
on which they connect is the root of the Higgs branch, and in order for it to touch 
the Coulomb branch classically, W'{— TOj) must be zero. A question is whether or not 
this last point might be modified in the quantum theory: while the metric on the 
Higgs branch does not get quantum corrected, this does not mean that the (position 
of the) boundary of moduli space can not get quantum corrections. In the study 
of the integrable system we usually take W'(— m) to be zero, however we have not 
found a special status of this extra condition on the superpotential: the Higgs branch 
continues to be there in the analysis. 

2.1.1 Quantum mechanics 

As is well known, the moduli space of supersymmetric gauge theory is governed by 
holomorphy and the non-renormalization theorems. The metric on the Higgs branch 
which was described above for instance does not get quantum corrected, see |13j . 
The quantum corrections to the Coulomb branch are quite well known since the work 
of Seiberg and Witten. Note that these arguments fail on the special points of the 
moduli space on which extra degrees of freedom become massless. Below, we will 
review the field theory analysis of the quantum corrections to the Coulomb branch 
found above. 
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2.1.2 Saddle point analysis 

First we want to study the theory in a regime where the theory is expanded around 
saddle points of the classical superpotential. We therefore want to fix the eigenvalues 
of the adjoint such that W vanishes for those values. Let us consider W'(x) = 
nr=i( x ~ a i)- ^ \ a i ~ a j\ >> 7^ j, then the saddle point analysis should give 
a good description of the theory. The gauge group is broken to U(Ni) subgroups 
centered at the a's, where each subgroup has an effective quadratic potential ~ /i$ 2 
for the adjoint scalar of this subgroup, and the masses of the hypermultiplets get 
shifted by a { . We can therefore reduce the analysis to the analysis of a Af = 1 theory 
perturbed by a quadratic potential with a certain number of massless fundamentals 
which can be combined into an analysis for the general superpotential. The vacuum 
structure now depends on the number of massless flavors charged under this subgroup, 
Mi, since this number determines the effective low energy degrees of freedom of the 
theory. These can be determined, for instance, by analysis of the anomaly consistency 
conditions (see e.g. [T%| HDj). We will analyze the system below mostly in the limit 
of vanishing masses and a large but finite /i. 

2.1.3 Mi < Ni 

The mass of the effective potential is very large (but finite), so the adjoint scalars can 
be integrated out of the superpotential (JIJ to give an effective theory for the meson 
field M\ = QiQ 3 a . Now add to this lagrangian the known N — 1 non-perturbative 
contribution of the Afneck-Dine-Seiberg |2Hj superpotential to arrive at the potential 

i Sfq-Mj _ x 

W = trM 2 + (Nt - M i )A i Ni ~ Mi det(M)^^7. (12) 

This potential is exact by the usual combination of holomorphy and symmetry argu- 
ments, supplemented by knowledge of weak coupling limits and analyticity, see e.g. 
[2*3] . The M = 1 scale A, can be related to the M = 2 scale by scale matching, 
]\3Ni-Mi _ ^N^2Ni-Mi_ ^y e can now ana ]y Ze the minima of the above potential by 

using the technique outlined in [23]. We find that the meson matrix can be diago- 
nalized by flavor rotations to a matrix with two different entries. Let the number of 
entries of one type, say ot\ be given by an integer r < [-£■]. In general we will have 
2N — Mi different solutions for every a±, leading to 2N — Mi different vacua at every 
r. If r = 4p then we have Ni — Mi/2 vacua. 

Bare quark masses can be added to (JT2J), which generically break the full flavour 
symmetry group. For a given r branch with massless quarks there are ( *) inequiv- 
alent ways of deforming this by giving masses to the quarks. The vacua obtained as 
minima of potentials of the form (j!2|) will be referred to as 'electric' vacua. 
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2.1.4 Mi = Ni 



In addition to the meson field degrees of (gauge invariant) freedom this case also has 
a baryonic degree of freedom. The effective M = 1 theory is described by a linear 
sigma model for the meson and baryons. Adding to this the potential inherited from 
the M = 2 theory we get an effective potential 

W = -^UrM 2 + K(det(M) — BB — A 27Vl ). (13) 

The vacua of this theory include the same type of vacua as for Ni < Mi when 
B = B = 0, so we can analytically continue the above counting results to this case. 
In addition to these, there is one vacuum with non-zero expectation value for the 
baryons B,B . 



2.1.5 Mi > Ni 

This case is fundamentally different from the previous one: the degrees of freedom in 
the low energy effective action are the Seiberg dual ^H] magnetic quarks and gluons 
for a gauge theory with dual gauge group U(Ni = Mi — Ni) and Mi flavors. The 
effective superpotential for these degrees of freedom is given by 

W d = qiM]q ] - ^trM 2 + [{Ni - M<)A • det(M)^^J , (14) 

where the part between parenthesis is only added when the rank of the meson matrix 
is Mj, since in this case the IR theory is a pure glue theory. This superpotential has 
two branches of solutions, depending on whether the meson matrix is degenerate or 
not (compare the classical analysis for the Higgs branches). If it is not degenerate, 
then the dual quarks can be integrated out and the analysis is the same as in the 
case Mi < iVj. In particular, we get the same number and degeneration pattern of 
vacua as we obtained there. We will refer to this case as 'electric' vacua. If the 
meson matrix is degenerate, the solutions are classified by an integer f, just as in the 
previous analysis, but which now appears for the dual quarks and mesons. We will 
refer to these vacua as 'magnetic'. The difference in the situation in the electric case 
is the coupling to the meson matrix. This coupling gives a mass to all quarks. In 
the IR we are therefore left with a classical gauge group U(N C — f) with no massless 
flavors which gives N c — f vacua with a surviving U(l). In addition to these two 
branches we also have a baryonic branch, just as in the electric case. 

The weak coupling analysis now yields the table ^ where we have only listed 
phases which appear in addition to the ones already found in the cases listed above 
it. The binomial factors are included to facilitate comparison to the results further in 
the paper. Note that the binomial factors are included and we stress that this is the 
result for zero effective quark mass. In this table we have abused notation slightly by 
labeling the baryonic branches by the integer r or f. 
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Table 1: Number of extra vacua for M = 2 softly broken to M = 1 by a quadratic 
superpotential. Here Aj = N c — Nf. 



case # vacua labeled by 

all ^ ( M r i ){2N i -M i ) r = 0,...,[M l }/2] 

(f)(N t -f) r = M t /2 

Mi > Ni 1 r = iV; 

M,>iVi (^)(JVi-r) r = 0,...,iVi 

1 r = A>- 



2.1.6 Strong coupling analysis 

The final step consists of analyzing the Coulomb branch dynamics encoded in the 
Seiberg-Witten curve of the gauge theory and especially on its singular points. We 
take the Seiberg-Witten differential X sw from 



X sw (z) = Z - (p\z) - \{P{z) - V(*))f^) dz = ^dln (y + P{z)) , (15) 



with 

.2 _ a \2N C -N f / , __, \N 



y l = P( z f - 4A 2Wc -^ (z + mi)"' =P 2 -AQ (16) 



the Seiberg-Witten curve. As noted in [21], there are several proposals for the exact 
form of the Seiberg-Witten curve. Specifically, there can be a term of form A 2Nc ~ N f 
times a polynomial of order Nf—N c if Nf > in the polynomial P. This corresponds 
to an inherent ambiguity in the definition of the quantum operators which reduce to 
-rtr$ fc in the classical limit for k > 2N C — Nf. In general there are several possibilities. 
Specifying a definition of these operators in the quantum theory fixes the ambiguous 
polynomial, see e.g. [22]- In this paper the natural definition by the gauge theory 
resolvent is used, 

Itr<$ fe >=i— / dzz k -\X sw , (17) 



k k 2ixi _ 

where the loop runs counterclockwise. Note that the Seiberg-Witten form also con- 
tains one power of z. 

Turning on a superpotential lifts most of the moduli space of the parent M = 2 
theory; the only points not lifted are those at which new degrees of freedom (dyons) 
become massless or which are fixed by non-zero expectation values of the quarks. In 
this paper we do not consider theories in which mutually non-local dyons become 
massless and we take all quarks to have the same mass m. Theories of the sort we are 
studying in this paper have been studied by various other methods in the literature, 
for instance by using geometric engineering 124.., Konishi anomaly techniques 



9 



exact potentials [2E], or by using a brane setup [2E]. The upshot is that the gauge 
theory vacua are described by a factorized Seiberg-Witten curve of the form 



P{x) 2 -AQ{x) = (x + m) 2r H 2 Nc _ r _ n+l (x)F 2n _ 2l (x) 
(W^\x) + f n+1 (x)) = B?(x)F 2n _ 2l (x), (18) 

where r labels the Higgs branches as before, if is a polynomial in x of order N c — 
r — n + I, B\ is a polynomial of degree / and W' n (x) is the derivative of the classical 
superpotential. 

In jT^j an extra contribution for n > 2N C — Nf to the non-singular part of the 
Seiberg-Witten curve was reported by revisiting the field theory analysis in j2Hj- The 
quantum operators used in that derivation however are not the quantum operators 
defined through the resolvent. In appendix El it is shown that if these quantum 
operators are defined through the gauge theory resolvent the extra contribution in 
the calculation disappears. 

The baryonic Higgs branch is much more easy, since the gauge theory is in this 
case fully Higgsed. The only way this can be reflected in the Seiberg-Witten curve is 
the case when the whole curve factorizes, y 2 = H 2 . 

In the following we will repeatedly use a partially integrated form of expression 
(|17|). Partial integration will give a boundary term in the loop integral because the 
logarithm has a branch cut which we choose to be on the negative real axis. However, 
this boundary term mostly cancels in calculations (as it should) since we will only 
calculate the pole at infinity. Here the function In (y + P(z)) ~ In (z N + l.o.) = 
\n(z N )+\n(l+0(l/z)). The last part is a nice holomorphic function in a neighborhood 
of z = oo. Now the boundary term cancels against the term coming from § z k ~ 1 In z N . 
The only exception to this canceling mechanism are terms 

±-jd{\og{z N )) = N. (19) 

These terms we will however not need in this paper since they correspond to the trace 
of the identity operator. We can therefore write 

~tr < $ fc >= ^ j dzz k ~ l In (J^{v(z) + P(z))j . (20) 
2.2 Compact ificat ions on R 3 x S 1 

Above we reviewed the vacuum structure of four- dimensional M = 1 gauge theories. 
Once we compactify on a circle, the vacuum structure remains the same, as explained 
in the introduction. However, the natural moduli of the theory are different. The 
adjoint scalar of the four- dimensional theory remains an complex adjoint scalar on 
M 3 x S 1 , but the 4d gauge field yields two additional scalars, one being the component 
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of the 4d gauge field in the direction of the S , the other being the dual of the 3d 
gauge field. Both scalars are compact. At a generic point on the moduli space of the 
M = 2 theory compactified on M. 3 x S 1 , the gauge group is broken to U(1) N , and 
the moduli space is therefore 4A-dimensional. In the analysis using the spin chain 
below, we will encounter complex variables Pi and qi. In the classical limit, they can 
be roughly thought of as the vacuum expectation values of the adjoint scalar $ and 
of the two additional scalars respectively. The precise relation of these variables to 
the quantum operators in the theory is something we will not discuss in detail, but 
will certainly suffer from an ambiguity similar to the one discussed in the previous 
subsection. In addition, from the integrable system point of view any set of variables 
p'i, q\ that is related to the original variables pi, qi by a canonical transformation is 
perfectly acceptable. We will see this explicitly in the next sections. Again, which of 
the canonical variables have the simplest interpretation in terms of the microscopical 
definition of the gauge theory is something which is outside the scope of the present 
paper. 



3 The classical SL(2, C) spin chain and its degen- 
erations 

In this section we review the classical inhomogeneous twisted periodic SL(2) spin 
chain. We follow mainly [27|. The curve with fundamental hyper mult iplets is a gen- 
eralization of the curve without, so the integrable system lurking behind should reduce 
in some limit to the periodic Toda chain. Since the spin-chain will be formulated in 
terms of 2 by 2 matrices, the periodic Toda system will first be written in terms of 
such matrices. 



3.1 Periodic Toda in terms of 2x2 matrices 



Periodic Toda (see e.g. [28] ) can be defined by an infinite-dimensional Lax matrix 
which commutes with the shift-by-N operator (is N-periodic). We can write this 
infinite dimensional matrix acting on an eigenvector as a recursion relation, 

AV"-^- 1 /^! + n / n + f n+1 = A/„, (21) 

which can be written as a matrix equation (in just 'site- local' variables for the matrix), 



fn+1 \ = (\-4>n -A 2 e-«" \ / f n 

e q -f n ) \ e qn J V e*-Vn-i 



(22) 
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Applying N matrices of this type will lead to the shift-by-N Toda transfer matrix 

Troda- 

Since our original infinite dimensional matrix was periodic with period N, we should 
focus on eigenvalues of the transfer matrix. These eigenvalues are given by the roots 
of the characteristic equation of this 2 by 2 matrix which is 

det(T Toda - w) = w 2 - tiT Toda w + det(T Toda ) = 0. (24) 

By defining y = 2w — trT^oda and using det {TtocLo) = A 2Ar we get the familiar form of 
the spectral curve, 

V 2 = (trT Toda ) 2 - AA 2N . (25) 

3.2 XXX 

The above formulation of the periodic Toda chain can be used to search for general- 
izations. Let us consider the following matrices (confusingly also called Lax matrices), 

L j (\) = \-\ j + y £s2 ) o*, (26) 

k=l 

T _ ( A - Xj + S{ S J 2 -iS J 3 \ . . 

3 V S( + iSi X- Xj - s{ J ' 1 ] 

where X{ are constants (trivial Casimirs) known as impurities and the are the Pauli 
matrices. We define a to be the identity matrix. For this system to be integrable 
the Lax matrices should satisfy quadratic r-matrix relations (see e.g. [2*§] . chapter 2), 

{L,(A) f Lj-(A')} = 6a [r(X - A'), L<(A) ® L 3 (X')] , (28) 

where we choose a rational r-matrix 2 (w-\) ELo cr a ®a a . This r-matrix is proportional 
to the interchange operator P : A® B — > B <8> A and the r-matrix relation is invariant 
under multiplication on the left or right of both Lax matrices with a constant matrix. 
With this r-matrix the Sj should satisfy a classical SL(2, C) Poisson algebra, 

{St,S b j } = i5 ab e ijk S a k . (29) 

Now define the transfer matrix and the associated spectral curve in analogy to the 
Toda case, 

i 

T(A) = J] Lj(\) (30) 

j=N 

det(T(A) - w) = w 2 - tiTw + det(T) = w 2 - P N (X)w + Q 2N (X) = 0, (31) 
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where the subscripts denote the degree of the polynomials. Defining y = 2w — trTxoda 
we obtain 

y 2 = (P N ) 2 - AQ 2N . (32) 

An important step in analyzing integrable systems is obtaining its conserved quan- 
tities. The spin chain variables have a natural quadratic Casimir, the spin length 
{S a ) 2 , which commutes with all S%. If we fix the values of the Casimirs, there 
remains a Poisson manifold of dimension 2N. The remaining variables can be ex- 
pressed in terms of N coordinates and momenta (see for instance appendix El) with 
the usual Poisson bracket. To completely solve the integrable system we therefore 
need determine N conserved quantities. Observe that products of Lax matrices at 
different sites obey the same r-matrix relations, e.g. {Li(A)L 2 (A) f Li(A')L 2 (A')} = 
[r(A — A'), Li(X)L 2 (X) <8> L 1 (\')L 2 (\')]. Taking the trace on both factors in the direct 
product now gives a conserved quantity of the spin chain. A convenient basis can be 
obtained by studying {T(A), T(A')}, which gives the coefficients of the trace of the 
transfer matrix as conserved quantities. 

Comparing with the Seiberg-Witten curve for U(N C ) with Nf = 2N C flavours 
in e.g. jTTj shows that the Q term should be (up to the modular functions of the 
conformal gauge coupling factor e"~, which we will discuss momentarily) Y\a=i (^+ m i)> 
where are the masses of the fundamental hypermultiplets 4 . In the XXX chain the 
polynomial Q 2 n is, 

N 

Q2N(\) = H((\-\i) 2 -K 2 ), (33) 

i=l 

where K 2 is the obvious quadratic Casimir of the Poisson algebra (|2~9j) at site i. We 
therefore obtain 2 masses from each site which are furthermore equal if is zero. 
Note that the full spectral curve of the spin chain only contains Casimirs of the algebra 
and is therefore itself a conserved quantity. 

As a final step in identifying gauge theory quantities, the modular functions have 
to be inserted in front of the Q, or for the more general (Nf < 2N C ) case introduce 
a scale into the problem. It has already been noted that the integrability condition 
(|28|) is invariant under multiplication of the Lax matrices by a constant matrix. We 
can therefore take generalized boundary conditions without affecting integrability, or 
in other words take a modified transfer matrix T = CT, where we take C to be 

In the case where Nf = 2N C , g is the modular function ^tjji ■ To really fit the curve in 

[T7j A should be shifted by gfi with fi = — ^- 5^^; and the masses should be redefined 

4 strictly speaking the rrii as used here are modular functions themselves. If we take a limit to 
the case of Nf < 2N C , the corresponding limit on the scale of the problem will reduce the m(r) to 
simple constant masses. 
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in the spin chain to be m = rrii — /i. However, since we will be reducing to non-finite 
gauge theories in a moment we will not do that here. In any case, in most cases 
considered in this paper the gauge theory scale will be introduced in another way 
than by the above C matrix, see subsection 13.3.11 

3.3 Degenerations 

In gauge theory the degeneration to lesser amounts of flavors is performed by inte- 
grating out massive degrees of freedom. On the level of the Seiberg-Witten curve this 
translates into taking several masses to infinity while keeping the product of these 
masses with the gauge theory effective scale fixed to a new effective scale. If our 
integrable models describe the gauge theory correctly, there should be a limit in the 
integrable system which reproduces this procedure. We will discuss how to take limits 
on the spin chain to integrate out masses. For one site and in different variables this 
was discussed in |2*7j . 

The SL(2, C) algebra at each site has a Casimir which we will want to keep fixed in 
this paper since it relates to the masses. Also the connection to the Toda chain is not 
clear in the variables Si, so we will re-express the SL(2, C) algebra in new variables 
q,p,K for each site. The Heisenberg-like variables q,p obey the standard Poisson 
algebra {qi,Pj} = <%■ An additional reason to switch to Heisenberg variables is that 
the 3d superpotential we will consider has known non-perturbative contributions in 
the pure glue sector of the form ^2 e ij+^-<n an d these appear naturally in Heisenberg 
like variables. Hence, these new variables are much closer to the naive gauge theory 
degrees of freedom. After some analysis (see the appendix (jXJ)) we arrive at the Lax 
matrix: 



where one of three possible forms is chosen, namely the one which is the most con- 
venient one for taking limits. The other two forms are connected to this one by a 
(possibly singular) canonical transformation. In this form the inhomogeneities Xj 
and the Poisson-algebra Casimirs Kj have been expressed in terms of the masses. We 
want to study degenerations to the Lax matrix of the Toda chain (|22|) which can be 
done by taking certain scaling limits. Immediately a problem presents itself: if we 
introduce the modular functions through a C matrix as above, then not all masses are 
multiplied by the coupling function in general which prevents us from taking masses 
to infinity directly. The origin of this problem is the site-locality of the masses at sites 
and the site global form of the coupling matrix. Since the masses are site-local, we 
will study a local limiting procedure below. However, this global versus local problem 
is the source of the existence of the different coordinate patches we will encounter 
later. 




(35) 
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3.3.1 Local limits 



As noted in the previous subsection, the integrable structure of the chain is invariant 
under multiplication by a constant matrix. Multiply (J33j) by constant matrices Aj = 

and Bj — ( ^ ) on the right and left respectively to obtain 



aj J " \ Pj y 

L= f X- Pj njm-,, • ,,,), \ 

V -l3j{m 2j -i + Pj)e q J (A + {m 2j + m 2j -i) + Vj)<*$j J ' 

At every site there are four different limits one can take: 

1. Take both otj and (3j to zero, keeping f3jm 2 j_i and ajm 2 j fixed (at A). We will 
call this limit (0) 

2. Take <x,- to zero and (3j to one, keeping a.jirn 2 j fixed (at A). We will call this 
limit (1"). 

3. Take (3j to zero and <x,- to one, keeping j3jm 2 j_i fixed (at A). We will call this 
limit (l d ). 

4. The trivial limit a = 1 and (3 = 1, which we will call (2) 

Note that the first limit reproduces the Toda chain 2 by 2 matrices. The gauge 
theory scale is introduced in the above by dimensional analysis and by comparing 
to the Toda chain. In principle we could also take an entirely different set of limits 
where in the end all the top left components of the Lax matrix would be put to zero 
when taking the Toda limit. However, it seems inconsistent to take such a limit on 
one of the Lax matrices and then take the above limit (0) on the next one: this 
procedure would for instance not reproduce a Seiberg-Witten curve of U(N C ) gauge 
theory without flavors in the case where all masses are integrated out. Since taking 
the above described limit or the alternative limit at all matrices will give the same 
spectral curve we will only consider the ones in the list above. 

We will use the following notation for limits: the total quark content of the chain 
will be given by a vector M := (e^, eN c -i, ■ ■ ■ , ei), where indicates the limit taken 
on a site. 



4 The proposal 

In this section we present our proposal for the exact superpotential of M = 2 U (N c ) 
gauge theory with Nf hypermultiplets in the fundamental representation on the man- 
ifold R 3 x S 1 , which is broken to M = 1 by adding a (polynomial) superpotential W 
for the adjoint superfield $. 
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4.1 Motivation 



As explained in the introduction, in order to write down the quantum superpoten- 
tial we merely have to find good holomorphic variables on the moduli space. The 
integrable system and its degeneration described in the previous subsection has all 
the required properties and therefore provides us with variables that are valid on 
at least an open subset of the moduli space. Therefore, our proposal simply boils 
down to using the conserved charges of the integrable system to build the quantum 
superpotential. There is a subtlety regarding exactly which conserved quantity we 
should associate to the operator tr($ fc ). As discussed in section 12.1.61 we will use to 
gauge theory resolvent to give a precise quantum meaning to tr($ fe ) and to define 
it unambiguously in terms of the conserved charges of the integrable system. This 
quantum definition is also the one for which the factorization of the Seiberg-Witten 
curve takes on the simple universal form of equation (18). 

4.2 The superpotential 

The form of the superpotential in 4d for the softly broken M = 2 theory is 

W = W(p) + v^QXQ? + V2Q^tQl. (37) 

In the pure glue case it was natural to replace the moduli space coordinates in the 
superpotential by (a specific basis of) the conserved quantities (Hamiltonian func- 
tions) in the integrable system, i.e. traces of powers of the Lax matrix. In the present 
case that would translate to replacing the moduli space coordinates by certain hamil- 
tonians of the XXX chain. In the direct application of this idea we then run into 
the problem of what to do with the Q's: since there are no dynamical variables in 
the chain which correspond naturally to the fermionic quark fields or the composite 
meson field variables. This is highly reminiscent of the lack of natural glueball field 
variables in the periodic Toda chain. Note that, just as the on-shell glueball fields can 
be calculated from the factorization of the curve, on-shell we can recover the meson 
vevs as well. Imagine integrating out the quarks, irrespective of whether this 
is allowed. We are then left with an effective superpotential of the same degree as 
the original one. However, in order to integrate out the quarks we need to choose a 
particular vacuum around which to expand. One would therefore expect to find only 
one branch, with fixed r, of the full vacuum structure. For an attempt to calculate 
an integrable system potential in a similar way see appendix O 

In the main body of this paper we will follow a more naive approach by taking 
iy($) and replacing the gauge theory Casimirs directly by a certain basis of the 
Hamiltonians of the spin chain. Perhaps not surprisingly, we will find that the equi- 
libria of the resulting superpotential generically describe several r branches at the 
same time. To make the proposal concrete we have to specify how to translate the 
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gauge theory traces in the superpotential to an expression in the conserved charges 
of the XXX chain. As stated above, and motivated by the strong coupling analysis in 
the gauge theory, we take the approach of using the gauge theory resolvent to define 
the operators Tr < $ fe > in the integrable system. This proposal reduces to the pure 
Toda chain one [5] in the Toda limit (on all the sites). 

In the following several vacua of this proposed superpotential will be studied. 
We will find that there can be different solutions to the equations of motion which 
yield the same Seiberg-Witten curve. These solutions therefore describe the same 
vacuum from the four dimensional point of view. However, on M? x S 1 they obviously 
correspond to inequivalent vacua. 

5 Addition and multiplication maps 

On field theory grounds some vacua of gauge theories with different matter content 
and/or different rank of the gauge group should be connected through the so-called 
addition and multiplication maps (see for instance JH]). In this section we show that 
these maps both have a natural counterpart in the integrable system. We will use 
that the e.o.m. derived from a polynomial superpotential can be written as, 



Here Oj is either pi or q^. Hence if the derivative of the transfer matrix with respect 
to the coordinates vanishes the equations of motion are satisfied. 

5.1 Addition map 

Suppose we have a solution of the spin chain e.o.m derived for a fixed superpotential 
W = ^2=0 ^tr($ fc ) = ^Ukgk in a N c ,Nf parent gauge theory. Using this it must 
then be shown that (some of the) solutions to the e.o.m of spin chain corresponding to 
a daughter gauge theory with the same superpotential and N c = N c + 1, Nf = Nf + 2 
can be generated. Likewise, there must be a solution in the N c , Nf which can be 
projected down to the parent theory. In the field theory analysis this map is 
known as the addition map and below we show how it has a natural interpretation 
in our proposal. 

Start with the transfer matrix of the Nf, N c theory and add one diagonal Lax 
matrix (x + m)l2 X 2 to it, which can be obtained by adding a spin chain Lax matrix 
with 2 masses fh and setting p — > —fh. Then it easily seen that the spectral curve 
of the resulting transfer matrix has the form (x + fh) 2 (P^ — 4Q) appropriate for a 
N c = N c + 1, Nf = Nf + 2 theory, with 2 flavours of mass fh on the Higgs branch. It 
then must be shown that the e.o.m. of the two theories have a simultaneous solution 




(38) 



17 



which has p — > — m, so that a solution in the parent theory can be lifted to the 
daughter. The generating function for the e.o.m. in the N c theory is obtained by 
adding a matrix at site N + 1 to a U (N c ) transfer matrix, 

trT- (*) = tr (( ,?~ *\ q {Th + P *>~ W ' ) ( hl h2 )), (39) 



^ T N C ( X ) = ( x ~ Pfijtu + (m + pfije q ^t 21 - 

(m + Pfi c )e q *ct 12 + (x + 2m + p^ c )t 2 2- (40) 

It is easy to see that the e.o.m. for o» G Pi,qi i ^ N c + 1 simply project down to the 
U(N C ) theory when pfj = —fh, since only the t's depend on these coordinates. In 
this case the e.o.m. can be derived from 

—=- — = {x + m)— — . (41) 

OOi OOi 

The equation of motion for q^ c is also solved by this value for p^ . The equation 
obtained by varying with respect to is quadratic in e 9 ^ on this solution since in 
that case, 

trT ~ t 22 - t n ~ e~ q ^t 2l - e q *ct 12 , (42) 

and the Seiberg-Witten curve is independent of e qfr <^ on the solution. The resulting 
equations can be solved for e 9 ^ which yields two solutions. The solutions for p^ 
and q^ c map two solutions to the e.o.m. of the daughter theory to the parent. These 
two solutions for q^^ correspond to the same Seiberg-Witten curve in the parent 
theory. Vice-versa, one solution to the parent theory e.o.m. can be lifted to two 
solutions in the daughter with equivalent Seiberg-Witten curve. It remains to be 
shown that the e.o.m. projected down to the parent theory can be derived from the 
same superpotential in both theories, but this follows trivially from (|38p. Note that 
on shell we have 

u k = -{m) k + u k . (43) 
5.2 Multiplication map 

The multiplication map can be derived in the 2x2 formalism just as it was for 
the pure glue gauge theory in the N c x iV c form in [S]. The derivation below is the 
extension of the integrable systems multiplication map to theories with fundamental 
matter. Consider a transfer matrix for a U(tN c ) gauge theory with tNf quarks which 
is the product of t transfer matrices Tjv c of the parent U(N C ) theory. 

It needs to be shown that solutions to the e.o.m. for the U(tN c ) can be projected 
solutions to the e.o.m. of the U(N C ) theory by periodically identifying coordinates 
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and momenta. The difference between the e.o.m. for qj and q a t+j can be studied in 
the U(tN c ) theory, with a an integer, 

r-r-)**-- < 44 > 

kj oq at+j J 

The transfer matrix can be written in a site-local form. Hence taking derivatives with 
respect to coordinates qj and q a t+j (or momenta) and periodically identifying with 
period t will yield the same transfer matrix in both cases. Therefore the right hand 
side of equation ()44j) vanishes and the e.o.m. of the U(tN c ) theory can be projected 
down to U (N c ) theory. 

It remains to be shown that the equations of motion can be derived from the 
same superpotential and how the Seiberg-Witten curves are related. Below we will 
first show that the curve in the U(tN) theory is related to the curve in the U(N) 
theory as 

Pw c ~ *QtN c = H 2 2(t _ 1)Nc (P 2 Nc (x) - 4det(T iVc )), (45) 

with H a polynomial of order 2(t — 1)N C . In particular, if one curve factorizes as 
in equation (fTHjl for a particular solution in that theory, so does the multiplication 
map of that curve when evaluated on the solutions which have been generated by the 
multiplication map of that solution. 

The trace of the product of transfer matrices, trT^, only depends on the two 
eigenvalues of T Nc , so the trace of the product of t of them can be written as a 
function of the trace and the determinant of T^ c . In fact, since for 2 by 2 matrices 

tr(T tNc ) = tr(T^) = tr^Jtr^ 1 ) - det(T Nc )tr(T N J- 2 ), (46) 

we have 

tr(Z*J 2tr(T^ e ) tr(T^) tr(T^ 2 ) 



(2det(r,v c ))*/2 2 v /det(T 7 vj2det(T 7Vc )(*- 1 )/ 2 2(det(T 7V J)(*- 2 V 
or defining C t = ^.^L and x = C x - v ' -- 



2 : 



(47) 



2(det(T Wc ))*/2 - - -1 2v / det(TjVc) 

C t (x)=2xC t _ 1 (x)-C t „ 2 (x), (48) 

which is the Chebyshev recursion relation (for Chebyshev polynomials of the first 
kind with the above initial conditions). Putting it all together we get, 

P tNc (x) = 2(det(T JVc ))'/ 2 C t ( ^L= ) , (49) 

\ 2 v /det(T J vJ J 
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as is expected on purely field theory grounds (see e.g. [30J) and which combines nicely 
with the determinant of T t N c to give the correct factorization of the curve. 



P&-4QtN. = 4det(T Nc y[c? 1 Pv ' U) 



2- v /det(T j 



= det(T N JY- 1 Ul 1 l-^^)(Pl(x)-4det(T Nc )). (50) 

Note that the first part of the last equation can be written as the square of a 
polynomial. This can be easily proven using the parity relations for the Chebyshev 
polynomials. Since the scale in the U(tN c ) theory must be real, det(T/v c ) must be a 
t-th root of unity and there are t inequivalent ways of lifting a solution of the U(N C ) 
theory to a solution of the U(tN c ) gauge theory. 

Using the above analysis of the Seiberg-Witten curve, it can now be shown that 
the equations of motion of parent and daughter theory in the multiplication map are 
derived from the same superpotential. From equation (J3%j) we get that 

PnM \ 6 Pi 



(~n I Pn(x) \ 5Pn_ 

XT = 9^7 f W & x (51) 

^^dct^J J V N ^ N 



The following relation can easily be proven using the generating functions for the 
Chebyshev polynomials, 

C' t (x)=tU t . 1 (x), (52) 

and it is then easy to see that the equations of motion in both parent and daughter 
theory are derived from the same superpotential. Note that there is no freedom to 
choose limits in these versions of the multiplication and addition maps, although when 
the above maps are combined with finite non-singular canonical transformations there 
certainly is such freedom. 



6 Examples 

This section contains several examples in which the vacua of several of our pro- 
posed superpotentials are calculated. Unless stated otherwise these examples are 
for a quadratic superpotential with W'(—m) = and with equal mass fundamental 
hypermultiplets. 
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6.1 N c = 2 
6.1.1 N f = 2 

We will do this case in some detail, to illustrate the (in)equivalence of the various 
limits that one can take on the full XXX chain. In the notation introduced in section 4, 
there are three different and inequivalent limits that one can take, namely (2,0), 
(1", 1") and (l u , l d ). Below we list for each of these three cases the value of P 2 (x), 
Q 2 {x), Ui = Tr($) and u 2 = ^Tr($ 2 ). The latter we compute using the gauge 
theory resolvent, so that it is slightly different from the naive quantity that one 
would extract from P 2 (x) (see equation (|138J) ). For simplicity we also introduce the 
notation y = e qi ~ Q2 . We find for the case (2, 0) 

P 2 (x) = (x -pi){x -p 2 ) - yA(p 1 + mi) - y~ l A{p 1 + m 2 ) 

Q 2 (x) = A 2 (x + m 1 )(x + m 2 ) 

ui = pi+ p 2 

u 2 = -ipl+p 2 2 )+A 2 + Ayip 1 + m 1 )+Ay- 1 (pi + m 2 ), (53) 
for the case (l u , 1") 

P 2 (x) = (x - pi)(x —p 2 ) - yA(pi + mi) - y^ 1 A(p 2 + m 2 ) + A 2 
Q 2 {x) = A 2 (x + m 1 )(x + m 2 ) 
Ui = pi+p 2 

u 2 = j (Pi + + Ay ( pi + mi ) + A y 1( yP^ + m 2), (54) 

and finally for (l n , l d ) 

P 2 (x) = (x - pi)(x - p 2 ) - yA 2 - y' 1 (pi + m 1 )(p 2 + m 2 ) + A 2 

Q 2 (x) = A 2 (x + nii)(x + m 2 ) 

ui = pi + p 2 

u 2 = -{pi +pj) + A 2 y + y' 1 (pi + m 1 )(p 2 + m 2 ). (55) 

As we expected, these three cases are related to each other by canonical transforma- 
tions. The second case is obtained from the first one by the transformation 

Pi^Pi-Ay~ l , P2^P2 + Ay^ (56) 

while the third case is obtained from the second one by the transformation 

y ^ Ay/(pi + mi). (57) 

Notice that this last transformation is not well-defined everywhere on the phase space. 
This illustrates the general phenomenon that different degenerations can in principle 
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describe different open subsets of moduli space, though on overlaps they are always 
related by a canonical transformation. 

The vacua of a purely quadratic superpotential (proportional to u 2 ) in the case 
of equal masses m\ = m 2 = m are easily found. In the third case (1", l d ), there are 
three types of solutions. The first are the r = solutions 

A i 771 

Pl =p 2 = ±A, y = — (58) 

for which 

P 2 (x) - 4Q 2 (x) = (x =f 2A) 2 (x 2 ± 4mA + 4A 2 ). (59) 

This follows the general factorization pattern summarized in (j!8|) . Secondly, we have 
the r = 1 solutions 



— m ± \Jm 2 — 4A 2 — m =)= \J m 2 — 4A 2 
Pi = g ' P2 = 2 ' V = 1 

for which 

P%(x) - 4Q 2 (x) = (x + m) 2 (x 2 - 4A 2 ). (61) 

This shows that both solutions give rise to the same Seiberg-Witten curve. They are 
nevertheless inequivalent solutions for the gauge theory on I 3 x S 1 . 

There is one more solution which is not strictly speaking a solution. It corresponds 
to the baryonic branch. Consider the e — > limit of 

Pi = P2 = ~ m + e, y — e/m. (62) 

One can easily check that the pi and g« equations of motion vanish in this limit, and 
that the Seiberg-Witten curve becomes 

P 2 2 {x) - AQ 2 {x) = (x + m + A) 2 {x + m - A) 2 . (63) 

However, this solution is at the boundary of the coordinate system. 

The r = and r = 1 solutions can be mapped directly to the first two limits (2, 0) 
and (l u ,l u ). However, the baryonic branch is more complicated, as this is exactly 
where the map between the (l u , l u ) and (l u , l d ) degenerates. One can take limits of 
the coordinates in (2, 0) and (l u , l u ) so that the Seiberg-Witten curve takes the form 
(jHSJ), but in none of those limits do the equations of motion also vanish. In any case, 
the baryonic root appears to be never a bona fide finite solution of the equations. 

6.1.2 Nf = 1 and 3 

For the limits (l u , 0) and (l d , 0) we find the expected 3 genus zero SW curves. For 
(2, l d ) we find 3 different phases (r = 0, 1 and a baryonic root). The r = 1 branch is 
triply degenerate, as expected. The baryonic root is located on the edge of the phase 
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space as above. The gauge theory analysis gives in this case 2 phases in addition to 
the three found. These are the phases obtained from the Seiberg dual quark potential 
in equation (fT4*|) when the meson matrix is degenerate. The non-degenerate meson 
matrix leads back to the electric vacua, which we do find here. 



6.1.3 N c = 2,N f = 4 

This is a simpler version of the case studied in [7] and for direct comparison we take 
the superpotential in this case to be purely quadratic W(x) = x 2 . The difference is 
that in that article the flavour masses were all assumed to be unequal and the solution 
was calculated in the weak coupling limit. We will not present the calculation here, 
but we find the expected 2 electric branches: two r = branches and a quadruply 
degenerate r = 1 branch. 



6.1.4 Unequal masses, Nf = 2 

We study this in the limit (l u ,l d ). Since this example is more involved we will 
elaborate slightly. We start by imposing a superpotential of the form W(x) = \x 2 . 
The e.o.m. for the momenta are easily solved. We are then left with one equation for 
the coordinates which, using the transformation e qi = bi, e Q2 = b 2 reads, 



b 2 (A 2 (V - b 2 2 ) 2 + h 2 (6i b 2 m 2 2 - (V + b 2 2 ) m 2 m 3 + b x b 2 m 3 2 ) J 



h 3 - h b 2 ' 



0. (64) 



This yields 5 solutions, a trivial one (b 2 = 0) and the four solutions of the quartic 
equation. The Seiberg- Witten curve for the trivial solution is 

y 2 Tiv = (x 2 - A 2 + m 2 m 3 + x (m 2 + m 3 )) 2 . (65) 

The solution to the quartic equation is not very insightful by itself and we were unable 
to show the factorization of the Seiberg- Witten curve directly. However, factorization 
can be shown through a different route 5 . Note that, since we are first and foremost 
interested in roots of the Seiberg- Witten curve, we can try to find it's double roots 
when evaluated on the solution for the momenta. We shall use the quadratic part of 
the numerator of equation (jMj) repeatedly to simplify expressions. Using this, it can 
be shown that the spectral curve simplifies to 

2 a a2 / , \ / i \ \ { 2 i bi(m 2 + m 3 )^ 2A 2 6i - b 2 



y 2 = -4 A 2 (x + m 2 ) (x + m 3 ) + [ x 2 + 1V , , *' x + ^ .(66) 

h + b 2 bi J 

We first try to find its double roots by studying the cubic equation %r = 0. Happily, 
this equation factorizes into a linear and a quadratic equation. The linear equation 



5 another route would be to compute the common divisors of the of the curve expressed in spin 
chain variables and the remaining equation of motion 
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has a root at 

(m 2 + m s )bi 

— (67) 

This root is a candidate for a double root of the spectral curve flfifij) . Plugging this root 
back into the curve we see that it is indeed a double root of that curve iff equation 
(JMJ) holds. We have therefore found the double root of the spectral curve. Using 
polynomial division gives the non-singular part, which has the simple form 

2 4A 2 fr 2 

(x) . (68) 

Numerical evaluation of the double roots shows that these consists of solutions which 
connect (at m 2 = m 3 ) to a doubly degenerate r = 1 branch and the expected two 
solutions for the r = branch. 



6.1.5 Cubic potential, Nf = 2 

We study a superpotential of the form W(x) = |x 3 — a 2 x, using the limit (1 M , l d ). 
Solving the equations for the coordinates yields three independent solutions. Two of 
these are related by A — ► —A. The third one yields two genus one Seiberg-Witten 
curves of the same form, 

y 2 = -4 (m + xf A 2 + (-a 2 + x 2 + 2A 2 ) 2 . (69) 

The extra contribution in the P-polynomial exactly cancels the leading term of the 
Q polynomial to form W'(x) 2 + / n _i just as expected from the field theory analysis 
in appendix El There are 2 main branches of solutions remaining with two remaining 
equations in each of them. Examination of these equations shows these have again 
3 sub-branches of solutions. Since the equations are related by A — > —A for the two 
main branches, we only have to explore one set of these 3 sub-branches. The first 
sub-branch gives again the above genus 1 curve (quadruply degenerate this time). 
The second sub-branch gives rise to 4 different genus zero curves and the last one 
consists of two doubly degenerate r = 1 curves. 



6.1.6 Cubic potential, Nf — 1 

We consider the same potential as above and the e.o.m. give in this case 2 main 
branches of solutions. One main branch consists of two sub-branches. One of them 
numerically corresponds to 6 genus zero curves, whereas the other sub-branch contains 
(numerically) a root at approximately x = —m. However, we were unable to obtain an 
analytical expression for these last curves. The remaining main branch corresponds 
to two genus one curves, as expected. 



24 



6.2 N c = 3 



6.2.1 N f =A 

We study the limit (l M ,2,l d ), since this case resembles the correct limit for N c = 
2, Nf = 2 and find 2 + 2 + 1 vacua for r = 0,1,2 respectively. These vacua belong to 
the electric phase. Furthermore we find the direct (addition map) lift of the N c = 2, 
Nf = 2 baryonic branch. It appears that we are missing the same (magnetic) vacua 
here as we did in the N c = 2, Nf = 3 case. The r = 2 branch yields 6 times the same 
curve, the two r = 1 branches four each, and the r = just has one solution for each 
curve, as expected from the field theory analysis. 

6.2.2 N f = 2 

We find 4 r = and 2 branches of the solution using the (2,0,0) limit. The r = 1 
branch (which has two solutions for each curve) corresponds to the addition map 
applied to N c = 2, Nf = 0, since it has the structure (x + m) 2 times the curve for 
N c = 2, Nf = with the same superpotential. The (l u , l d , 0) limit shows one extra 
phase for a solution of the e.o.m. which is on the edge of the phase space (qi — > oo). 
This fase has Seiberg-Witten curve factorization 

y 2 = (x + m)\(x + m) 2 + 2A 2 )((x + mf - 2A 2 ). (70) 

which is a genus 1 curve. This is an example of a solution which is not expected from 
the gauge theory analysis which lies on the boundary of the phase space. 

6.2.3 Nf = 3 

For this interesting case we checked several different limits. Just as in the N c = Nf = 2 
case we generically find the expected 3 r = and 3 (with 3 solutions each) r = 1 
electric branches. Only limits which contain a link like in the Nf = 2 = N c case, i.e. 
(1", l d ) also produce the baryonic root. 

6.3 N c = 4 

In order to check the multiplication and the addition map we study the four flavoured 
theory. The multiplication map can be studied by taking the limit (l u , 1^, 1", l d ). This 
gives 4 phases at r = 0, 4 at r = 1, and 2 phases at r = 2. Actually, the phase at 
r = 2 has a (x + m) 6 piece in the Seiberg-Witten curve, but if we shift the mass 
in the superpotential from + mx to a slightly different mass m we see that only 
4 roots of the Seiberg-Witten curve remain at x — —m. In addition we find the 
multiplication map lift of the baryonic branch. However, in addition to these phases 
we also find a Seiberg-Witten curve of the form ((x + m) 8 + A 8 ) and one of the form 
(m + y) 4 ((y + m) 2 — A 2 )((y + m) 2 + A 2 ) which are genus 1 curves, both very similar 
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to equation (JTHj) . Just as in that case, these solutions are on the edge of the phase 
space. In fact, the genus 1 curves seems to be the addition map lift of this 'wrong 
curve'. We checked also a different limit for this case, which yields the same phases 
with the exception of the baryonic root. 



6.4 Nf < N c 

We take the limit (Q Nc ~ N f(l u ) N f), N c to be larger than to 2 and the superpotential 
to be W(x) = |x 2 + mx. We take the masses to be all unequal. We compute the 
superpotential, 

W = -J2p* + J2 Aeqi ~ 1 ~ 9l ^ + m ^ + Yl A2 e ?w ^ + m^ ft . (71) 

i=l i=l i=Nf+l i=l 

To see this is indeed the correct expression for the quadratic potential note the fol- 
lowing: for Lax matrices at a site which contains masses, the constant lower right 
matrix element only contributes to the trace of the transfer matrix at order x Nc ~ 3 . 
Since a quadratic potential is needed that can be calculated from the coefficient of 
x Nc ~ 2 and x Nc ~ x that element can be safely ignored for calculational purposes. The 
general form of the quadratic superpotential can then be inferred by comparing to 
the Toda case. From equation (fTTj) we easily obtain and solve the equations of motion 
for p, 

Pi + m = -Ae qi - 1 - qi i = l,...,N f (72) 
Pi + m = i — Nf + 1,...,N C (73) 

Her the variable go is defined as q^. The equations of motion for q are a bit more 
involved, 



e 9i-i-9i(p. + m .) = e « + i(p. +1 +m i+1 ) i = l,...,iV7 - 1 

e qN f- l ~ qN f{ PNf +m Nf ) = Ae QN f~ qN f +1 i = N f 

Ae 9 *- 1 -" 1 = Ae qi ~ qi + 1 i = N f + 1 . . . N c - 1 

\ e lN c -i-g Nc = e qN c -qi( pi+mi ) i = N c , 



(74) 



or using the solution for P and the definition e 94_1 qi — Zi, 

Ziirtii - Azi - m) = z i+1 (m m - Az i+1 - m) i = l,...,N f -l 
ZNf(nT<Nf — AzN f —m) = Az^ f+ i i = Nf 

Zi = z i+ x i = N f + 1 . . . N c - 1 

Azn c = z\(ra\ — Az\ — m) i = N c . 

(75) 

The variables Z{ obey the constraint Yi z i — 1- At this point it is important to note 
that the equations of motion indeed have a solution, up to the constraint. Also, in this 
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case all solutions give finite coordinates. Evaluated on the solutions for the momenta 
the transfer matrix becomes 



Nf+1 / x , m \ 1 / x , _, , m y ( m i y m \\ 

T (x) = n ( x ) n ( 1 i T ~ zj ~^> y ^ 

i=N c ^ ^ ^ ' j=N f ^ ^ Z 3 J 

Here we have obtained the variables Zj by inserting '1' between Lax matrices 

as 

1 ° V 1 ° ) (77) 
and realizing we are only interested in the eigenvalues of the transfer matrix. 



6.4.1 Equal masses 

First take all masses equal to m. Then Nf simple quadratic equations have to be 
solved, which are solved for Nf less than N c as = €iy/—ZN c , with q equal to either 
+ 1 or —1. The constraint now reads, 



i N ' {z N f^ 



N f 



n^=i. (78) 



i=i 



For even Nf this constraint can be solved straightforwardly for zjy c for any choice 

of £j. We therefore have a total number of 2 N f(N c — Nf/2) solutions. When Nf 

i 

is odd, we have to square the constraint which gives z^ c = (— 1) 2iVc_JV / . Inserting 
this back into the constraint then fixes one of the e, arriving at the same number of 
solutions. On the basis of the field theory analysis we expect f jv^O (N c — Nf/2) + 

ES'o 2 " 1 ("') (2iV c - N f ) solutions for even N f and £fi /2 ~ 1/2 ( N /) (2N C - N f ) for 
odd Nf, which can easily be shown to be equal to the number of found solutions 
using 2 k = £i=o (i)- ^ Nf equals N c the above analysis can be repeated by solving 
everything in terms of Zi = | ^tv c | (note that there are N — 1 epsilons). Then we see 
that there are 2N C solutions by squaring the constraint, and the constraint fixes one 
of the epsilons for the two possibilities (zn c ) Nc = ±1. For the rest of this subsection 
we set the mass m for convenience to zero and consider Nf < N c since the case 
Nf = N c is a simple generalization of the argument below. The mass can be restored 
by shifting x — > x + m and shifting p — > p — m in the solution. 

The above counting argument gives the total number of solutions for all r branches. 
Below we show how to obtain the factorization of the curve. Use the solution of the 
Zj in terms of the z^ and (cyclically) redistribute factors of J~zn c in (fTBj) to get, 

Nf+l / T \ 1 / x 



T(x) = A N TT A V . TT S lV c ■ V c • (79) 



N c * ' j=N f 
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Define x = , p — , so 



^Z+i / \ 1 



m = ( ^A)»tr n ( * ;)n( s tr i)- < 8o » 

i=N c v 7 j=JV> v J 7 

This expression for the transfer matrix consists of building blocks 

-Ui) H*-"^) < 8i > 

These matrices share a basis of eigenvectors and commute. In fact, 

A + ■ A- = xB (82) 

This immediately leads to the identification of the r-branches in this solution: we 
take a set of Nf — r A + 's and r A~ J s, r < Nf, zn the appropriate solution of (J7H|) 
and apply ()82)1 to get the transfer matrix 

T(x) = (^zV c A) N (x) r B Nc - N f +r (A ± ) N f +r } (83) 

which will yield a spectral curve of the form y 2 ~ x 2r G(a;). The ± depends on whether 
r < [-^] or not. From equation (f?8*J) we see that we have N c — ^ solutions of this type 
for even Nf and r < ^f. Therefore for r < ^£ there are 2 (^) (iV c — ^) solutions 
for each r. The factor of two comes from the fact that for fixed r both the solutions 
with r v4 + 's as the solutions with r v4~ contribute. For r = there are therefore 
only N c — solutions. For odd Nf we have to take into account that there has to be 
either one A + or A~ more than the other. Then we have f 1 ) (2N C — Nf) solutions 
at each r-level, so the counting of states works out level by level. 

Now writing the building block matrices in an eigenvalue basis gives a transfer 
matrix with diagonal elements: 

N c -N f / , \ Nf—r 



Tn(x) = (J^FA) N I | + \/ (|) 2 - 1 



2 V 2 




f-* + V ( f )2 - 1 I (S1) 



Nf-r 



T 2 2{x) = (V^TAY 



2 




55) 



The square of the difference between the diagonal elements is the discriminant of the 
characteristic equation of the 2 by 2 matrix, so is the Seiberg-Witten curve, 

y 2 = (T n -T 22 ) 2 . (86) 
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By first reducing by using (}82|) it is easily seen that the remaining Seiberg-Witten 
curve has the form polynomial times \J (|) 2 — 1 and thus factorizes. 

Motivated by we show how to write this in term of linear combinations of 
Chebyshev polynomials. Chebyshev polynomials of first and second kind obey the 
identities 

Tn(x) = \((x+ vV - 1))" +{x- vV - 1))") (87) 
U n ^(x)Vx^l = vV ~ 1))" vV - 1))") • (88) 

It is easily seen by expanding 

(| + . ± ' = ±Q (| ± v^i) ' • («» 

that the trace of the transfer matrix is a linear combination of Chebyshev polynomials 
of the first kind, while the Seiberg-Witten curve is a linear combination of Chebyshev 
polynomials of the second kind times the expected non-singular piece. 

6.4.2 Unequal masses 

In this case far less specific information than above can be obtained since to obtain 
more concrete results in this case the constraint equation F] — 1 has to be solved, 
which is hard to do in general. However, note that inserting the solution for the Zi in 
term of z^ c into equation (|76|) leads to a transfer matrix of the form 

N f+ 1 / \ i 

T(x) = A N J] ( * + % ZN o ) n ( X+^+f **> ) , (90 ) 

i=N c ^ ' j=N f ^ Zj ' 

where zj are the roots of 

Zj ( A ~ Zj ~~ 1^ = ZNc ' ^ 
Set x = -p7= and £• = to obtain 

^-^"n(* ;)n(*!*i). 

The building block matrices in this transfer matrix all commute so can be diagonalized 
simultaneously, as before. The eigenvalues now differ from previous cases only by the 
value of a in equation (fKHJ) since the differences in the Lax matrices between the cases 
with and without equal masses are proportional to the identity matrix. Hence we 
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can (again) write the spectral curves of these transfer matrices in terms of Chebyshev 
polynomials and the spectral curve factorizes in the appropriate manner, 




From the point of view of gauge theory it is interesting that we find a concrete 
realization of the polynomials studied in (20] generating the factorization for both 
equal and unequal masses for the full range of Nf < N c . 



6.5 Massive vacua 

In line with the calculation in j2] we show how to generalize the above computation 
to the massive vacua of an arbitrary (single trace) polynomial superpotential for 
Nf < N c . Since the factorization into a genus zero curve shown above is very special 
the solution of the integrable system equations of motion shown is a natural starting 
point to look for a solution for all polynomial superpotentials W = Ylt-^i 9k\- From 
the gauge theory point of view for massive vacua the Seiberg-Witten curve must 
factorize into a genus zero curve which we take to be 




According to equation (fT%|) there must then be polynomials G and / such that 

Y 1 V'"; ,(•''! K(rf ■ (9--,) 




holds. Since this factorization problem is exactly the same as in the pure glue case 
which was analyzed in \5\ we only state here that for such polynomials to exist, the 
constant c must be such that the coefficient c$ in 

1 - 

W{y/WM£ + 7) - c) = J2 c ^ ( 96 ) 

i=—n 

is zero. We will first show below that a slight generalization of that solution indeed 
solves the equations of motion up to one equation which we then show to be equal to 
the equation expected from the gauge theory. 

We take our solution of the equations of motion in (jHEJ) to be equations (J73j) and 
(J75|) together with the constraint Yi z i = 1> where in equation (|73p we replace m by 
the variable c. We have already shown that solutions of this type yield factorized 
Seiberg-Witten curves with non-singular part as in equation (}95|) . but we still have 
to show that these solutions actually solve the equations of motion. 
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First change variables as x = y/z^hx — c in the integral of equations ()38j) above. 
The equations of motion for the g« are satisfied because the appropriate derivative of 
the trace of the transfer matrix evaluated on the equation of motion gives (after some 
manipulations) 

ST N r — A ( -v^vT 





T, 



N-l; 



(97) 



5qi " ~ v -\/ z n c 

where T/v-i is the transfer matrix obtained by eliminating Lax matrix Lj, evaluated 
on the solution to the equation of motion. In particular, Tjv-i is diagonalized a known 
basis of eigenvectors. In this basis the other matrix becomes 





-\fZN~A 



^A 








\{x - \lx 2 - 4) 







(98) 



Therefore the trace of equation (|9*7jl vanishes and the e.o.m. for the q's are satisfied. 
The e.o.m. for the momenta can be derived from the transfer matrix 



ST, 



N 



5pi 



-1 Zi 





T) 



N-l- 



(99) 



Changing basis to the basis of eigenvectors of the other Lax matrices in the transfer 
matrix (and some manipulations) yields 



-1 Z{ 





Vx 2 -4: 

Zi+\{-x+^x 2 -4) -z l +\\-x-\[x 2 ^i) 



z l +\{x-y/x 2 -4) 



(100) 



Vx 2 -4 Vx 2 -4: 

Therefore the trace of equation (}99|) becomes 

ST N (2 Zl + x){1$_ x - T%_ x ) - y/W=Z{T$_ x + 7f_ x ) 



tr- 



V. 



x 2 



with 



T, 



N-l 



mil 
A N-l 





m22 

1 N-l 



(101) 



(102) 



Now we can also express the Seiberg-Witten curve on the solution in terms of the 
transfer matrix Tjv-i using the fact that on that solution 



Ti 



N 



\/x 2 — 4: 







Zj + \(x+yJx 2 -A) 
Vx 2 -4 



rpll 

1 N-l 







rp22 

1 N-l 



(103) 



As explained above, the Seiberg-Witten curve is the difference between the diagonal 
values, which gives something proportional to equation (jlOlj) . We therefore arrive at 



5P_ 

Spi 



1 



solution 



(104) 



31 



Hence the equations of motion for all the momenta collapse into 1 equation. Com- 
bining this computation with the shift in x and changing variables x — ► £ + | yields 

j> W'(y/zFM£ + j) - c) j# = 0, (105) 

which is equivalent to equation (J5fi|) expected from the gauge theory point of view. 
We have therefore proven that this particular solution solves the equations of motion 
for the integrable system and is such that equation (fTHjl holds. Note that this holds 
for all Nf < N c and for arbitrary masses. 

As a technical aside, note that in the Toda case the derivatives of the transfer 
matrix with respect to the momenta can be summed such that on the solution the 
following equation holds: 

^dP _ dP 
^r 1 dpi dx' 

This reproduces, when evaluated in the resolvent, exactly the condition in that 
trW'(M) should vanish, where M is the periodic Toda chain Lax operator. 

It is instructive to realize at this point that the factorized Seiberg-Witten curves 
obtained in this subsection should be the full field theory answer to the factorization 
problem for massive vacua of gauge theories with fundamental matter with Nf < N c . 
In particular, they account for all possible vacua in a unified way. However, the 
above construction is constructive in the sense that we do not prove that these are all 
possible factorizations. On the other hand, they are all possibilities that are expected 
from the weak coupling point of view by a simple counting argument. To compare 
to other approaches to gauge theory we calculate in the final subsection the vevs for 
the adjoint scalar in these vacua. 



6.6 Adjoint scalar vevs for massive vacua 

The vacuum expectation values for the adjoint scalar in the massive vacua discussed 
above are given by 



<u k >= ~ cf> a*" 1 tog ( VJ w - V — I dx. (107) 



For the massive vacua 



P{x) + ^P^AQ = 2T U = 2(VSj£A)"« J] ( | + z, + ^ - 1 J , (108 
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with x = Jlitf^ . The expectation values can be split into a part which also appears 
in the pure glue case and extra terms depending on the flavours, 



<C Uk > <C Ilk ^pure glue ^ Uk ^matter 

1 



<C U}~ ^pure glue 



<C Itfc ^matter 



2vri 



2vri 



X^A^lOg 



a; 



V 



x fc 1 log I 1 + 



(109) 



rfx (110) 



§+v(§) - 1 . 

Calculating the glue term first we get, by splitting the logarithm into 

1 



log ( 1 + - ] + log | - + 

X, 



1 



(112) 



2 V 4 x 2 

two finite integrals which can be performed. The first vev is easily calculated to be 

<^l>purcgluc= ~CN C . (113) 

The needed contour integral of both logarithms in equation ()112|) can be evaluated, 
with the help of the expansion 6 



1_ f2l\ l_ 



i=i 



to give 



\k] 

Nc^fk\ ( 

<Mfc>pureglue- fc ^ I % J I 

1=0 v 7 v 



2\U „\k-2l 



(114) 



;ii5) 



The matter contribution can be calculated by using the change of variables x — > £ + ^- 
As an intermediate result 



1 

'2iri 



(£+7^(1-4)108(1 + 



£=oo 



E 

m=0 



1 



j — m + 1 lm 



1 f_ y .V'-2m+l 



(116) 



is obtained by straightforward computation. Plugging this result into the needed 
integral we arrive at 



<C Ufo ^matter 



N f fc-i [i] 

i=l j=0 m=0 



fc - 1 



\j-2m+l 



my j — m + 1 



(117) 



6 Takc a derivative on the left hand side w.r.t. y = \ and observe that this is — J[^\ y with C(y) 
the generating function for the Catalan numbers 
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The full answer for the adjoint scalar vev in the massive vacua is now the sum of 
equations (|115jl and (|117jl . Comparing above expressions to the ones obtained by a 
matrix model calculation in [SZj we immediately identify 



c 



-T 



R. 



(118) 



Their constraint is exactly the constraint Yl z i = 1- This can be seen by multiplying 
the two different solutions to equation (jHH) for each site with a mass. For the solutions 
of that equation a consistent sign has to be chosen to connect to the matrix model 
result. With these identifications the relation to equations (67), (68), (69) and (70) 
in j3Zj is obvious. The relation with the remaining equality (71) is more obscure, but 
for k = 2 there is perfect agreement. It would certainly be interesting to check the 
matrix model calculation further against the above field theory result. 

6.7 Toward a general proof of factorization 

In this final subsection we indicate how a proof of factorization of the Seiberg-Witten 
curve can be obtained up to a technical detail. In previous work jH] we relied on a 
theorem in [2E] which relates straight line flows on the Jacobian of the spectral curve 
to the equations of motion of the periodic Toda chain. This connection is lacking in 
the case of the classical spin chain, although the fact that it is integrable and periodic 
strongly suggests that the reasoning should hold in exactly the way. In this paper 
however, we are mainly interested in the critical points of the potential as can be 
obtained by using equation (JHBJ), reproduced here for convenience, 



Assuming that the transformation Jacobian is not trivial, the equations of motion 
amount to the vanishing of certain one forms on the Riemann- surface, a condition 
actually identical to the conditions derived in [6 . Therefore, modulo the technical 
assumption, this amounts to a proof of factorization as in equation (fl~8|) of the spectral 
curve of the integrable system in exactly the same way as derived from the gauge 
theory. 

7 Integrating in the mesons 

A natural question in this framework is whether or not one can integrate in the 
mesons. A natural reason to suspect that this might be the case is the fact that the 




These equations can be written as 




(119) 
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solutions we found for Nf < N c are calculationally very similar to vacua of the po- 
tential (|12|1. The gauge theory interpretation of the solutions of both superpotentials 
is also similar. In addition, pure field theory arguments indicate that it might be 
possible to integrate in the mesons and even the baryons [33]. Finally, in [3H] it was 
shown that the minima of the ADS potential and the minimum of the matrix model 
potential agree for a quadratic potential. In this section we show how to directly 
integrate in the meson fields for a quadratic potential in the case Nf < N c and how 
this reproduces both four and three dimensional results. 

The quadratic XXX superpotential (obtained from equation (|7T|). reproduced here 
for convenience) is for the limit 1" on the first Nf sites 

/ N c N f 

w=n \ j2pi+ 2 J2 Ae9l+1 ~ 9l ( pi + 

\i=l i=l 

Here qN+i is identified with q\. A convenient way to rewrite this is in terms of new 
variables y\ = Ae qi+1 ~ qi (pi +m»). These satisfy one constraint which we impose using 
a Lagrange multiplier S, and we get: 

W = /i$>, 2 + 2y 4 ) + Slog ( >± ) . (120) 




i=i 



In this equation A only appears in the logarithm and it couples to S. Therefore, 
S is to be identified with the familiar glueball field of pure Af = 1 theories. We 
now introduce new Lagrange multipliers Mj for the masses and new fields rrtj and we 
rewrite ()120|) as 

Nc fA 2N c~ N f V\ Nf (v +m )\ Nf 
W = + 2 ^) + Shg[ >± + Y^irm - mm (121) 

i=i V iiVj J j =1 

If we integrate out M we regain ()120|) . On the other hand, we can also integrate out 
Pi, rhi and all y except yN c and we find 

Nf ( M 2 \ f N 2Nc ~ N f2 Nc LL Nc \ 

W/ = E(-^ + «M) + ^ HN.-N t -l)S + S log ( gWjSW ^J ■ 

(122) 

Integrating out S gives 

Nf f M 2 \ / \2N c -N f nN c „Nc\ N c -N f -i 

W = E [~t + ^«»MN<-N f -l) ( yw-nJf ; ) ■ d23) 

From now on we identify in the above expression FJ Mj with the determinant of the 
meson matrix. In order to compare to the three-dimensional M = 2 theories studied 
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in jHSj, we want to take the limit \i — > oo, keeping A 2Nc ~ Nf /j, c fixed by scale matching. 
Furthermore we set the quark masses to zero and arrive at the superpotential for the 
four dimensional gauge theory we study, compactified on a circle 

W = Y + (N c - N f - 1) (Y det M) ^-"f- 1 , (124) 

with a convenient normalization. If Nf — N c — 1, then integrating out S in the final 
step yields the right quantum constrained moduli space. 

Another direction from equation (j!22|) is integrating out i/n c and S, which yields 
equation (|T2*j) upon putting the masses to zero and trivially redefining \i as 2/x. In 
this case we can also treat Nf = N c and arrive at the correct quantum constrained 
moduli space. Baryons can be integrated in by using the procedure described in one 
of the examples in [3.3], but this does not seem to yield illuminating results. 

When N c < Nf < 2N C in principle one can also apply equation (jl21j) by redefining 
the y variables appropriately. However, it is difficult to obtain a concrete general 
expression for the superpotential in this case. For the low order cases we checked (up 
to N c = 3) the answer is exactly the potential expected for the electric branches of 
the theory. We hope to return to this issue elsewhere. 

8 Discussion 

We have shown in this paper that the XXX spin chain provides good holomorphic 
coordinates for part of the moduli space of U(N) gauge theory with Nf fundamentals 
compactified on M. 3 x S 1 . Some coordinate patches touch the baryonic root at the edge 
of phase space, but there are sometimes other solutions to the equations of motion 
on this edge as well. These other solutions do not have a gauge theory interpretation. 
Being on the edge neither class is a bona fide vacuum of our superpotential, but it 
is certainly suggestive and it would be interesting to understand this further. Since 
the coordinates provide a coordinate patch for part of the moduli space it is natural 
to wonder whether there is a more general integrable system which provides coordi- 
nate patches for the other parts of moduli space, or even coordinates that are valid 
everywhere. In particular it would be interesting to find a system which includes the 
baryonic root holomorphically. It might be worthwhile to see if there is a more gen- 
eral formulation of the integrable system starting from the formulation in SL(2, C) 
variables. 

We have also shown that it is natural to integrate in the mesons into the superpo- 
tential by Legendre transforming with respect to the masses. The baryonic operators 
cannot be naturally integrated in and it remains an open question whether or not 
there is an extension of the integrable system in which there is a parameter multi- 
plying the baryon. This parameter would then correspond to the variable b in [HH] 
and could be used to integrate in the baryon. More generally it remains a question 
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whether the meson and dyon vevs have an interpretation in the integrable system. 
An avenue of approach to this problem starts with noticing that the meson vevs are 
simply the value of the matrix model resolvent at x — —mi. This gives at least on 
shell a method of backtracking to the integrable system what these vevs are. 

Yet another open problem remains the microscopic interpretation of the phase 
space variables pi and qi in the three dimensional gauge theory. A way to examine 
this issue could be to study perturbative instanton contributions in the gauge theory 
to find the coordinate transformation of our results to the right microscopic variables, 
just as was used in jS]. 

As indicated before, we expect that the factorization of the Seiberg-Witten curve 
on the minima of our superpotentials which is expected on field theory grounds can 
be proven along the lines of the second proof in [Sj. We have indicated above how 
such a proof works in principle, but it would be interesting to close the small technical 
gap in that proof. Interesting generalizations of results in this paper include but are 
not limited to studying other gauge groups, for which the integrable system is, to our 
knowledge, unknown and/or to higher dimensional gauge theories. 
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A Heisenberg variables 

In this appendix we will rewrite m term of the Heisenberg variables for which 
{Qi,Pi} = 1- The aim will be to obtain something which resembles equation (|22jl . We 
will drop the indices on the g's, p's and S"s to declutter notation. First rewrite the 
algebra as, 

{S ,S ± } = ±S ± {S + ,S_} = 2S (125) 

We will begin with the Ansatz Sq = — p — A, and read (|125|) as a system of differential 
(derivational) equations. The first can be easily solved by taking, 

S± = f±(p)e ± « (126) 

where the functions f± contain a constant of integration. Putting this into the second 
equation we obtain, 

U(p)e- q {e«, /_(p)} + e«/-(p){/+(p), e~*} = 2(-p - A,) (127) 
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Take the functions f± to be polynomials of degree 2, f± = a± + b±p + c±p 2 . This 
yields three different solutions, 



2Ai^ _ 1 ^2 

/+ = a+ 

/_ = a_ + b-p 

U = (-2t + a t)~l 

f- = a- 



f- = aJ^n -p 



U = a + -^p-±p 2 



(128) 
(129) 
(130) 



The constants can be connected to Aj and Ki through S 2 + S + S_ = Kf, which yields, 

f = X ?~ K ? 2 A, 1 2 

f+ = a+ 

/_ = 6_(±^ i + A,) + 6_p 

/+ — V b- ) ~ b-P 

f- — a ~ 

f K}-\} 2 \i 1 2 ( 133 ) 

where in the second line 2 solutions are possible (which differ only in which mass is 
in which off-diagonal entry). 



B Seiberg-Witten curve factorization in the strong 
coupling approach 

In this appendix it is shown how to calculate factorization loci of the Seiberg-Witten 
curve when turning on a classical superpotential in the gauge theory approach, par- 
ticularly when the superpotential is of high degree (W ~ tr$ n+1 , n > 2N C — Nf). In 
this appendix single trace superpotentials are studied for Nf < 2N C . 

The theory at the root of the r-th Higgs branch [13] is &U(r) x U(l) Nc ~ r gauge 
theory, with Nf massless hypermultiplets charged under the U(l) subgroups of U(r). 
On special points along the root l—r extra massless hypermultiplet degrees of freedom 
appear which are charged under the U(l) subgroups. At these points the Seiberg- 
Witten curve factorizes as 

l-r 2N-2l-2r 

P(x) 2 -4Q(x) = (x-m) 2r l[(x-p t ) 2 J] (x-q j )=H l (x) 2 F 2N _ 2l (x). (134) 

i=i j=i 

Here P(x) = E£o s ^ c ~* and Q = ^ 2Nc ~ Nf (e£o rh t x N f-^ . 
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We will denote the extra massless dyons (and their scalar components) by di and 
dj. Note that since r eigenvalues of $ are fixed we are (classically) left with a U (N c —r) 
gauge group. For this theory we can write the superpotential 

l— r n+1 

W = V2^2<i i a i di + ^ f g k u k . (135) 

i=l k=l 

Since the theory is classically (for a generic point on the Higgs branch part of the 
vacuum) broken to a U (N c — r) gauge theory with Nf — 2r flavours, the u k correspond 
to multi-trace operators for k > N c — r. Therefore, for n < N c — r the analysis below 
works 7 . 

There is an ambiguity in the definition of the operators u k , since for higher values 
of k there can be several quantum operators with the same classical limit. We will 
denote the 'classical' u k by u k in the following. In other words, u k = | Xlj a i f° r 
P(x) = Yli( x ~ a i)- However, for the u k appearing in the superpotential ()135|) we will 
use the quantum operators whose vacuum expectation values are obtained from the 
gauge theory resolvent. The quantum vevs u k can be expressed as 

The first logarithm gives for k < N c exactly the Newton-Girard relations between 
power sums and the Sj, which can be written as symmetric polynomials of the eigen- 
values. The second term produces the quantum corrections. Clearly, in the limit 
A — ► 0, Q — > and vfjp —>■ vfj^ as expected. We define u k = u k + Au k - 

Before proceeding to the general case, we first calculate the first few corrections. 
Since only residues at v — are needed, the quantum contributions Au k from the 
second logarithm can be expanded as 

a«, - — Idw-^ («™*-»t f!^M + \^ Nf ( v^om y 

AUk ~ 2<Ki X [ V (v^P(l/v)y + 2 V \(v»'P(l/v))*J + 

+ o(v 6N °- m f)J. (137) 

The fraction in the integrand can now be Taylor-expanded around v — since it does 
not contain poles any more. The quantum u k and the classical u k differ only when 

7 For n — N c — r we do not have to add dyons in the case the Seiberg-Witten curve does not 
degenerate further. In this case the superpotential can be minimized by writing it first in terms of 
the N c — r eigenvalues. Then they can be taken to be equal to the N — r different roots of W such 
that the non-singular part of the curve is W' N _ r (x) 2 + Jn } -2t- For n < 2N C — Nf this has the 
expected form of W'(x) + f n -i- A more refined argument based on the gauge theory resolvent is 
also possible. 
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k > 2N C — Nf and the first two corrections are 

k = 2N C - N f Au k = A 2Nc - N ? (138) 

k = 2N C - N f + 1 Au k = A 2Nc ~ Nf (m x + 2*i) - ~A. mc ~ 2Nf 5w c -N f ,i- (139) 

z 

The point of this calculation is that it matters which superpotential is used in (|135|) 
in a calculation along the lines of j2.6j to calculate Seiberg-Witten curve factorization 
loci. We will do this below for a (single trace) superpotential derived from the gauge 
theory resolvent. Equation ()135j) may be varied with respect to u k to obtain these. 
However, in that calculation we would need — %r. This is hard to calculate since the 



or 

,(<?) 



Sj are highly non-trivial functions of the u k . It is easier to vary equation ()135|) with 
respect to the classical moduli, since part of the calculation is then standard. This 
yields 

n =V2E dA^h- (140) 



Repeating the steps in |26j gives, 



*M = B 'm +0 & (14i; 



with 



n+l n+1 n+1 r. (q) 

w\x) ^y.^- 1 ( 142 ) 

r=l 1=0 r=l " U r 

B q is a polynomial 8 of degree less than or equal to I — 1 and to get the neat first term 
in W'(x) we have used the assumption that the degree of the potential is less than 
or equal to N c . By comparing leading powers in x in equation (j!41j) . it follows that 
N c — n<l<N c — l. Squaring equation ()141|) now gives an equation akin to equation 

(3.44) in ma, 

W>l(x) = BlF 2n . 2q + f n ^(x) + ^fjT 1 + OX). (143) 

n 2l\ X > X 

This equation indicates that B q F 2n -2q — ~2N C +N f with g a polynomial of 

the indicated order. From other considerations we expect to find in the end W'^ (x) 2 + 
/ n _i and below we show that this is actually the case if the extra contributions on 
both left and right hand side of this equation are taken into account. 

The left hand side of the equation can be rewritten in terms of the classical 
superpotential and its quantum contributions using equation ()142|) . Using the obvious 
definition W' = W'n + W'^ the left hand side of the equation reads 

(W' {c) ) 2 + 2W' {c) W' iq) (x) + {W l{q \x)f. (144) 



3 in [23 this polynomial was denoted by H(v) 
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Note that the extra terms on the right hand side appear at the same order as those 
on the left hand side. The extra terms on the left can be calculated by starting with 

where P denotes ^jjr- The last term in the integrand can be expressed as, 



-Nc-j 



which gives after some massaging, 



(146) 



5 v^JVc -j 

The second term above is suppressed in the integrand by a factor of x~ Nc ~ l and will 
therefore not contribute to the residue integral as long as the order of the superpo- 
tential n does not exceed N c . Plugging this result back into W'™(x) yields 



n+l n+l 

W'^Hx) = 



EE^- 1 ^/, ^-i= — (148) 

1=0 r=l Zm Jx'=oo Jl_ 4Q +1 _4Q \ ^ J 

For functions of the form f(x) = YliLo a i xP ~ l 

f( x )\ + = E T~ / a^VOO**' (149) 



2?ri 

i=0 



X '= oo 



holds where f(x)\ , = Y^i=o aixP % an d ^ ne ' res idua at infinity integral' is taken to be 
the same as in equation (fTYj) . Defining r — 1 = n — j, 



n+l n „ 4Q 

W' {!) (i) = Ef'E lH rf -p 2 da / ( 150 ) 



gives 

/ 1 /4Q 



W'^faO = W' {c) \ 



x 



1 4Q _|_ i 4Q V -P 2 

- 1 - p2 ~T -L p2 



(151) 



Therefore we arrive at an expression for the left hand side of equation (|14Hjl 

4Q 

W' {c \x) 2 + 2W'} c) W ,( - q \x) + {W ,( - q \x)) 2 = W ,{ - C \xf + [ W' {c \x^ 2 1 



(152) 
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The field theory moduli coordinates on the left hand side of equation (|143J) are the 
double roots, while on the right hand side these are the Sj. These are related by the 
factorization of the Seiberg-Witten curve and using this the extra term on the right 
hand side of equation (j!43|) is 



4Q(x)B 2 



(W^(xY + g m ). (153) 



In the above equation m denotes the order of the extra polynomial in the fac- 
torization formula. To start of it is of order 2n — 2N C + Nf, which is smaller than 
2n. This polynomial order was derived from equation (j!43|) . However, we see that 
the leading term which does not depend on the polynomial cancels on both sides, so 
the polynomial is of one order lower. Then the argument can be repeated until g is 
of order n — 1. This completes the argument that even for high order superpoten- 
tials the factorization of the Seiberg-Witten curve is captured in equation (fT%j) if one 
consistently uses the gauge theory resolvent to calculate the quantum operators u^. 



C An observation 

A standard technique in matrix model calculations for theories with fundamental 
matter is integrating out the matter content: the resulting matrix model can then be 
tackled by the usual techniques. The modified potential takes a form like 

~ log(det($ + m)) (154) 

for every integrated flavour. Considering that the translation for the pure glue case 
proceeded by directly identifying the matrix integral variable with the Lax matrix of 
the periodic Toda chain, we might try to postulate that one has to take the modified 
matrix model potential and insert the Toda chain Lax matrix into it to obtain the 
effective action on IR 3 x S 1 . This would be in line with the observation in [3~3] that 
matrix model calculated superpotentials for systems with matter in the fundamental 
can be constructed from data of systems without. Note that in this way we can 
introduce at most Nf < N c flavours. Also, this kind of potential leads to equations 
of motion of high order rapidly. We calculated (some of) the minima of a potential 
W = mx + \x 2 + h log(det(a; + m)) for N c equal to 2 and 3 to check this idea. Note 
that we leave a term h of dimension 2 in front of the potential undetermined for 
reasons which will become clear below. In the matrix model this is Nf. 



N c = 2 

We plug the periodic Toda chain Lax matrix into the potential. In a full calculation 
we would have to use the field theory resolvent to calculate the potential correctly, 
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but here we will follow a more direct route: we simply take the ^-independent part of 
the determinant before taking the logarithm. Solving the equations of motion we get 
6 curves. Four of these have a double root. In addition, there is a doubly degenerate 
curve of the form 

y 2 = ((x + mf + 2A 2 + hA 2 )((x + m) 2 - 2A 2 + hA 2 ) (155) 

Now if this curve is to have double points, h has to be ±2, where it has a double root 
at —m. The rest curve has the form of (x + m) 2 ± 4A 2 . Note that the calculation 
with the XXX chain gives a minus sign in this curve. The remaining 4 curves can 
now be evaluated for the same two values of h, and these give 3 different curves since 
2 solutions collapse to the same curve solution: 

y 2 = ((i + m ) 2 T 4A 2 )(a; + m) 2 (156) 
y 2 = (x + m)(m + x + 4:VTlA)(m + x + 2V^[A) 2 (157) 
y 2 = ( x + m)(m + x -4v / =FlA)(m + x - 2 v / =FlA) 2 (158) 

Note that we do not get the correct form of the curve for the two solution which 
probably should correspond to the r = branch. Also, we get one more curve at 
r = 1 than is expected from the field theory analysis. To further check this type of 
construction we turned to N c = 3. 

X c = 3 

We first have to fix the coefficient h. For this purpose we use the solution with equal 
momenta p± = P2 — Pz- This gives twelve solutions to the equations of motion, all 
of which have 2 double roots in the Seiberg-Witten curve (for unfixed h even). This 
branch has a curve with double roots at — m if and only if h = 0. We therefore try to 
find a h such that it reproduces at least one of the double roots of the r = 1 branch 
of the 3 flavoured theory. We find that we can fix h such that 2 of the curves have the 
correct double roots (but not the single roots), but then we can numerically simply 
verify that the other curves do not have double roots at a recognizable location. This 
indicates that this approach to finding the factorization locus of Seiberg-Witten curves 
with matter in the fundamental in this particular setup does not work. It might be 
that the solutions of the equations of motion we have not considered do have the 
correct double roots. However, then one would still have to explain that there is a 
huge excess of solutions, for which there is no (intrinsic) reason to discard them. We 
have therefore checked that the constructed potential for three colors cannot (fully) 
describe a system with 2 flavours, as one might suspect. 

Discussion 

Although an intriguing observation at the level of two colors, the constructed potential 
in this particular form obviously does not reproduce four dimensional field theory 
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answers for theories with more than 2 colors and is therefore not the correct integrable 
systems-inspired exact superpotential for systems compactified on a circle. However, 
based this observation we conjecture that there is a (generic for all N) non-polynomial 
superpotential for the periodic Toda chain which reproduces Seiberg-Witten curve 
factorizations for supersymmetric gauge theories with fundamental matter (Nf = N c ) 
which for 2 colors resembles the superpotential displayed above. Using the field theory 
resolvent to calculate the non-polynomial potential might be a good starting point 
to improve the above calculation. In any case, a further study of equilibria of non- 
polynomial potentials for the Toda chain would certainly be very interesting. 
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